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Abstract. We describe a powerful method for determining the equation of state of 
an ultracold gas from in situ images. The method provides a measurement of the 
local pressure of an harmonically trapped gas and we give several applications to Bose 
and Fermi gases. We obtain the grand-canonical equation of state of a spin-balanced 
Fermi gas with resonant interactions as a function of temperature [T| . We compare our 
equation of state with an equation of state measured by the Tokyo group in [2 , that 
reveals a significant difference in the high-temperature regime. The normal phase, at 
low temperature, is well described by a Landau Fermi liquid model, and we observe a 
clear thermodynamic signature of the superfiuid transition. In a second part we apply 
the same procedure to Bose gases. From a single image of a quasi ideal Bose gas we 
determine the equation of state from the classical to the condensed regime. Finally the 
method is applied to a Bose gas in a 3D optical lattice in the Mott insulator regime. 
Our equation of state directly reveals the Mott insulator behavior and is suited to 
investigate finite-temperature effects. 
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1. Introduction 

Ultracold gases are a privileged tool for the simulation in the laboratory of model 
Hamiltonians relevant in the fields of condensed matter, astrophysics or nuclear physics 
[3]. As an example, thanks to the short-range character of interactions, ultracold Fermi 
mixtures prepared around a Feshbach resonance mimic the behavior of neutron matter in 
the outer crust of neutron stars [HIS]. For cold atoms, the density inhomogeneity induced 
by the trapping potential has long made the connection between the Hamiltonian of a 
homogeneous system and an ultracold gas indirect. Early experimental thermodynamic 
studies have provided global quantities averaged over the whole trapped gas, such as 
total energy and entropy |6i i7j, collective mode frequencies [8], or radii of the different 
phases that may be observed in an imbalanced Fermi gas [9l [101 E]- Reconstructing 
the equation of state of the homogeneous gas then requires to deconvolve the effect of 
the trapping potential, a delicate procedure that has not been done so far. However, 
the gas can often be considered as locally homogeneous (local density approximation 
(LDA)), and careful analysis of in situ density profiles can directly provide the equation 
of state of a homogeneous gas [12], [13], [1], [Hj. In the case of two-dimensional gases, 
in situ images taken along the direction of tight confinement obviously give access to 
the surface density |T5l [161 [171 [18] and thus to the equation of state |T9]. For three- 
dimensional gases, imaging leads to an unavoidable integration along the line of sight. As 
a consequence, inferring local quantities is not straightforward. Local density profiles can 
be computed from a cloud image using an inverse Abel transform for radially symmetric 
traps [2^. A more powerful method was suggested in [13] and implemented in [HHl]: as 
explained below, for a harmonically trapped gas the local pressure is simply proportional 
to the integrated in situ absorption profile. Using this method, the low-temperature 
superfluid equation of state for balanced and imbalanced Fermi gases have been studied 
as a function of interaction strength [H [E]. In this paper we describe in more detail the 
procedure used to determine the equation of state of a spin-unpolarized Fermi gas in the 
unitary limit [1]. We compare our data with recent results from the Tokyo group [2], 
and reveal a significant discrepancy in the high-temperature regime. In a second part 
we apply the method to ultracold Bose gases. From an in situ image of ''Li, we obtain 
the equation of state of a weakly-interacting Bose gas. Finally, analyzing experimental 
profiles of a Bose gas in a deep optical lattice [21], we observe clear thermodynamic 
signatures of the Mott insulator phases. 

2. Measurement of the local pressure inside a trapped gas 

In the grand-canonical ensemble, all thermodynamic quantities of a macroscopic system 
can be derived from the equation of state P = f{fi,T) relating the pressure P to the 
chemical potential fi and temperature T. P can be straigthforwardly deduced from 
integrated in situ images. 

Consider first a single-species ultracold gas, held in a cylindrically symmetric 
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harmonic trap whose frequencies are labeled oj^ = ujy = uj^ in the transverse direction, 
and ujz in the axial direction. Provided that the local density approximation is satisfied, 
the gas pressure along the z axis is given by |13] : 

n/^.,T) = ^n(z), (1) 

where n{z) = J dx dy n{x,y, z) is the doubly-integrated density profile, Hz = jJ^^ — 
\mu:1z^ is the local chemical potential on the z axis, /i° is the global chemical potential. 
n{z) is obtained from an in situ image taken along by integrating the optical 
density along x (see Fig{T]). As described below, if one independently determines the 
temperature T and chemical potential /x", then each pixel row of the absorption image 
at a given position z, provides an experimental data point for the grand-canonical 
equation of state P{^z,T) of the homogeneous gas. The large amount of data obtained 
from several images allows one to perform an efficient averaging, leading to a low-noise 
equation of state. 

This formula is also valid in the case of a two-component Fermi gas with equal spin 
populations if n{z) is the total integrated density. The method can be generahzed to 
multicomponent Bose and Fermi gases, as first demonstrated on spin-imbalanced Fermi 
gases in [T|, [H] . 

3. Thermodynamics of a Fermi gas with resonant interactions 

In this section we describe the procedure used in [1] to determine the grand-canonical 
equation of state of a homogeneous and unpolarized Fermi gas with resonant interactions 
(a = oo). We also compare our data with recent measurements from the Tokyo group 
P, [2]. We then study the physical content of the equation of state at low temperature. 

3.1. Grand- canonical equation of state 

In the grand-canonical ensemble, the equation of state of a spin-unpolarized Fermi gas 
in the unitary limit, can be written as 

P(/x,T) = pW(^,T)/^T(C), (2) 

where P*^"-* (/x, T) is the pressure of a non-interacting two-component Fermi gas and 
C = exp{fx/kBT) is the inverse fugacity. Since P^'^\^,T) is known, the function /it(C) 
completely determines the equation of state P(/i, T). Let us now describe the procedure 
used to measure it. The pressure profile of the trapped gas along the z axis is directly 
obtained from its in situ image using equation ([T]). One still has to know the value of the 
temperature T and global chemical potential fi^ in order to infer hxiC)- We use a small 
amount of ^Li atoms, at thermal equilibrium with the ^Li component, as a thermometer. 
We then extract /i° from the pressure profile, by comparison in the cloud's wings with 
a reference equation of state. For high-temperature clouds {ksT > we choose /i*^ so 
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Figure 1. Scheme of the local pressure measurement: the absorption of a probe beam 
propagating along the y direction provides a 2D image on the CCD camera. Integration 
of this image along x provides the doubly-integrated density profile n{z) and, using 
equation ([Ij, the pressure profile along the z axis. 

that the wings of the pressure profile match the second-order virial expansion [22] (see 
FigJ2}i): 

For colder clouds, the signal-to-noise ratio is not good enough, in the region where ([3]) is 
valid, to extract /i° using the same procedure. We thus rather use the equation of state 
determined from all images previously treated as a reference, since it is accurate on a 
wider parameter range than ([3]) (see FigJJb). We then iterate this procedure at lower 
and lower temperatures, eventually below the superfluid transition. By gathering the 
data from all images and statistical averaging, we obtain a low-noise equation of state 
in the range 0.02 < C < 5 (see FiglSK). 

3.2. Canonical equation of state 



In [2] a canonical equation of state E{n,T) expressing the energy E as a function of 
density and temperature was measured using fits of absorption images taken after a 
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Figure 2. Determination of we plot the data from an in situ image as P/2kBTX'^g 
versus —fi/kBT = V{z)/kBT — fi^/kBT (black points). A wrong choice of /i° 
corresponds in this representation to a translation of the data in abscissa. We adjust 
/i° so that the wings of the pressure profile match a reference equation of state (in red), 
(a) For high-tcmperature clouds, we use the second-order virial expansion Q (b) For a 
lower temperature pressure profile, we minimize its distance with the averaged equation 
of state deduced from higher temperature images (in red) in the overlap region. 

short time-of-flight. In situ density profiles were deduced by assuming a hydrodynamic 
expansion. The temperature was extracted from the cloud's total potential energy at 
unitarity, using the experimental calibration made in [7]. In FigjSb the data from [2] 
is plotted as E{n,T) / E^'^\n,T) as a function of 9 = T/Tp, where n is the total atom 
density, Tp is the Fermi temperature, and E^'^^ (n, T) is the energy of a non-interacting 
Fermi mixture. 




Figure 3. (a) Grand-canonical equation of state of a two-component Fermi gas with 
resonant interactions from [1 (black dots), (b) Canonical equation of state from the 
Tokyo group 2 (open circles) and from the ENS group (black dots) . The dashed black 
line is the ideal gas equation of state, the dot-dashed (solid) black line is the second- 
(third-) order virial expansion, the solid green line is the Fermi liquid equation ^ 
and the solid blue line is the fit function ^ in the superfluid phase. The superfluid 
transition occurs at C = 0.05. 
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The comparison between the two equations of state requires to express our data in 
the canonical ensemble. The density n = dP/dulT is calculated by taking a discrete 
derivative, and we obtain the black points in FigJSb- While the two sets of data are in 
satisfactory agreement in the low-temperature regime T/Tp < 0.4, they clearly differ in 
the high-temperature regime. The disagreement of the data from [2] with the second- 
and third-order virial expansions calculated in [221 123] indicates a systematic error in this 
regime. This is possibly due to a breakdown of hydrodynamics during the time-of-flight 
as expected at high temperature. 

3.3. Fermi liquid behavior in the normal phase 

Above the superfluid transition and in the low-temperature regime 0.05 < C < 0-5, our 
data is well modelled by a Fermi liquid equation of state 

where C,n = 0.51(1) and m* = 1.12(3)m respectively characterize the compressibility of 
the normal phase extrapolated to zero temperature and the effective mass of the low- 
lying excitations. The agreement with is better than 5% in a large parameter range 
0.33 /i < ksT < 2 fi. Our value of ^„ is in agreement with the variational Fixed- 
Node Monte-Carlo calculations C,n = 0.54 in [2l], ^„ = 0.56 in [25j, and with the 
Quantum Monte-Carlo calculation ^„ = 0.52 in [26]. It is surprising that the quasi- 
particle mass m* is quite close to the free fermion mass, despite the strongly-interacting 
regime. Note also that this mass is close to the effective mass m* = 1.20 m of a single 
spin-down atom immersed into a Fermi sea of spin-up particles (the Fermi polaron) 
[27ll25l[28l[29l|30l[T2l[IIlE]. 

3.4- Superfluid transition 

The deviation of the experimental data from (j4]) for ( < 0.05 signals the superfluid phase 
transition. This transition belongs to the U{1) universality class, and the critical region 
is expected to be wide [31] in the unitary limit. Assuming that our low-temperature 
data belongs to the critical region, we fit our data with a function 

P(/i, T) = P^L(/., T) + A{Q - O'-'^HiCc - C), (5) 

where H is the Heaviside function and a ^ —0.013 is the specific heat critical exponent, 
measured with a very good accuracy on liquid ^He ^2]. We obtain the position of 
the superfluid transition (c = 0.05, or kBTc/ fi = 0.33, in agreement with the value 
kBTc/fi = 0.32(3) extracted in jl] using a simpler fit function. We thus confirm more 
rigorously our previous determination of the superfluid transition. In the appendix 
we discuss the validity of local density approximation around the superfluid transition. 
Under our current experimental conditions, the deviation from LDA is very small. 
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4. Thermodynamics of a weakly-interacting Bose gas 

In this section we apply equation ^ to the case of trapped Bose gases. First we test the 
method by determining the equation of state of a weakly-interacting Bose gas [33], [M] . 
We use an in situ absorption image of a ^Li gas taken from [35] (see FigJlK). ^Li atoms 
are polarized in the internal state \F = l,mF = —1), and held in an loffe-Prit chard 
magnetic trap with u)r/2TT = 4970 Hz and uJz/2tt = 83 Hz, in a bias field — 2 G. 
Thermometry is provided by a gas of ^Li atoms, prepared in |-F = \,m.F = and in 
thermal equilibrium with the ''Li cloud. 




z (mm) 



Figure 4. (a) Integrated density profiles n{z) for tlie ''Li component (blacks dots) and 
the ^Li component (open circles). The solid line is a fit of the ®Li component with a 
finite-temperature Thomas-Fermi profile, yielding T ~ 1-6(1) /iK. (b) Thermodynamic 
function g{C,) determined from the ^Li profile. The solid line is a fit of the data with 
a Bose function in the non-condensed region and a mean-field equation of state in the 
condensed region (see text). The dashed line is the equation of state of a classical 
gas g{C,) = C^^. The difference between the dashed and solid lines around C = 1 is a 
consequence of Bose statistics. 



4.1. Determination of the equation of state 

The equation of state of a weakly-interacting Bose gas can be expressed, in the grand- 
canonical ensemble, as: 



where ( = e~^^'^^^ is the inverse fugacity and XdsiT) = ^J2'nh? /rnksT is the thermal 
de Broglie wavelength. The pressure profile is calculated using ([T]). We aim here at 
measuring g{C)- We obtain the global chemical potential value = 0.10 /c^T by fitting 
the ^Li profile in the non-condensed region \z\ > 50 /im with a Bose function: 
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Combining the measurement of the pressure profile, of the cloud's temperature T and 
global chemical potential we obtain the thermodynamic function g{C) plotted in 
FigSb. 

4-2. Analysis of the equation of state 

In the region C > 1 the data agrees with the Bose function g{Q = (75/2 (C) expected for 
a weakly-interacting Bose gas. The departure from the thermodynamic function of a 
classical gas g{() = and especially the fact that g{() > 1 above the condensation 
threshold, is a thermodynamic signature of a bosonic bunching effect, as observed in 
[36| l37t 138] . The sudden and fast increase of our data for ^ < 1 indicates the Bose- 
Einstein condensation threshold. In the local density approximation framework, the 
chemical potential of a weakly-interacting Bose-Einstein condensate reads: 



where my is the ''Li atom mass and 077 is the scattering length describing s-wave 
interactions between ''Li atoms. We neglect here thermal excitations in the condensed 
region. Integrating Gibbs-Duhem relation at fixed temperature dP = nd/i between the 
condensation threshold Q and ( < (c, and imposing the continuity at C = Cc we obtain 
the equation of state in the condensed phase: 

9(0 = 9^/2{Q + ^^(log^ C - log' Cc). (6) 

Fitting our data with the function g{Q given by (jS]) for C < Cc and with 5'5/2(C) fo^' C > Co 
we obtain Cc = 1-0(1) and 077 = 8(4)ao = 0.4(2) nm. The uncertainties take into account 
the fit uncertainty and the uncertainty related to the temperature determination. The 
condensation threshold is in agreement with the value Cc = 1 expected for an ideal Bose 
gas, the mean-field correction being on the order of 1% [391 SO]- Our measurement of 
the scattering length is in agreement with the most recent calculations 077 = 7(1) Oq 

Extending this type of measurement to larger interaction strengths on Bose gases 
prepared close to a Feshbach resonance would reveal more complex beyond-mean-field 
phenomena, provided thermal equilibrium is reached for strong enough interactions. 

5. Mott-insulator behavior of a Bose gas in a deep optical lattice 

Here we extend our grand-canonical analysis to the case of a *^Rb gas in an optical 
lattice in the Mott insulator regime. By comparing experimental data with advanced 
Monte Carlo techniques, it has been shown that in many circumstances the local density 
approximation is satisfied in such a system |12]. We analyze the integrated density 
profiles of the Munich group. Fig. 2 of [21j . 
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5.1. Realization of the Bose-Huhhard model with ultracold gases 

Atoms are held a trap consisting of the sum of a harmonic potential Vh{x,y,z) and a 
periodic potential 



created by three orthogonal standing waves of red-detuned laser light at the wavelength 
A = 27r/k = 843 nm. The atoms occupy the lowest Bloch band and realize the Bose- 
Hubbard model ^3] : 



with a local chemical potential /i(r) = — Vh{r). The index i refers to a potential 
well at position r^, J is the tunneling amplitude between nearest neighbors, and U is 
the on-site interaction, U and J being a function of the lattice depth [5]. The slow 
variation of Vh{r) compared with the lattice period A/2 justifies the use of local density 
approximation. 

We consider here the case of a large lattice depth Vq = 22Er, for which J ~ 
0.003 ?7 ~ 0, and assume that the temperature is much smaller than U. In this regime 
the gas is expected to form a Mott insulator: in the interval G [(p — l)U,pU], where p 
is an integer, the atom number per site remains equal to p, and the density is equal to 
n = p{2/X)^. Integrating Gibbs-Duhem relation between and /i, we obtain that the 
pressure P is a piecewise linear function of /i: 



5.2. Determination of the equation of state 

We use a series of three images from |2T], labeled a, b and c, with different atom numbers 
= 1.0 X 10^ Nb = 2.0 X 10^ and = 3.5 x 10^ (see FigJSji). The integrated profiles 
n{z) are not obtained using in situ absorption imaging but rather using a tomographic 
technique, providing a ~ 1 fj.m resolution. The pressure profile is then obtained using 
equation ([T]). 

Each image i = a,b,c plotted as P as a function of — |ma;^z^ provides the equation 
of state P(/i) translated by the unknown global chemical potential By imposing that 
all images correspond to the same equation of state (in the overlapping fi/U region), we 
deduce the chemical potential differences between the different images /i^ — /i^ = 0.56 U 
and — = 0.61 U (see FigJSb). Gathering the data from all images, we thus obtain 
a single equation of state, translated by which is still unknown. We fit this data with 
a function translated by from the following function, capturing the Mott insulator 
physics: 



Vq (sin^(/cx) + sin^ (ky) + sin^ {kz 



))^ 




(7) 




where (p — 1)U < fj, < pU. 



P 



for /i < 



f/(A/2)-3 
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Figure 5. (a) Integrated density profiles n{z) corresponding to image a (open squares) 
b (black dots) and c (crosses), from [5T]. (b) Determination of the global chemical 
potential difference /xj? — by superposing the equations of states given by each image. 



rzi-^ for < ;U < Sfii 

6lll II — 5 III r r r e- 

ni-^ + n2 — — — tor d^i < ^ < o^i + 0^2 

ni— + n2— + n3 tor diii + 6^2 < fi, 



with 5/i2, ni, 77-2 and as free parameters. The value = 1.51 U yielded by 

the fit thus corresponds to the condition P ^ when ^ 0. Once it is determined, we 
obtain the equation of state of the Bose-Hubbard model in the Mott regime, plotted in 
FigH 



«o=0 



ny=0.9(l) n2=2.0(l) "5=3.1(1) 




Figure 6. Equation of state of a Bose gas in an optical lattice, in the Mott insulator 
regime. The solid line is a fit with a piecewise linear function capturing the Mott 
insulator behavior. The slope dP/d/i provides the density in each of the Mott zone, 
m = 0.9(1), n2 = 2.0(1), 713 = 3.1(1). 
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5.3. Observation of a Mott-insulator behavior 

After fitting the value of the other parameters resulting from the fit exhibit the 
characteristic features of incompressible Mott phases. The occupation number in the 
first Mott region is Ui = 0.9(1) atom per site and the size is 6fii = 0.9(l)f/. The second 
Mott region occupation number is n2 = 2.0(1) and its size is 6fi2 = l.l(l)f/. Finally, 
the third Mott region occupation number is = 3.1(1). These values agree with the 
theoretical values Ui = i and = U , in the T = and J = limits. 

5.4. Estimation of finite temperature effects 

The equation of state deduced from the experimental data is also suited for investigating 
finite-temperature effects. Since sites are decoupled in the regime J <^ f/, ksT 
considered in this study, the finite-temperature equation of state is easily calculated 
from the thermodynamics of a single site [HI H5] : 



Fitting now the experimental data with ^ and T and /i° as free parameters, we deduce: 



This value is in agreement with a direct fit of the density profiles and number statistics 
measurements [30] • This temperature is significantly smaller than the temperature 
ksT* ~ 0.2 f/ at which the Mott insulator is expected to melt |13]. Second, this 
temperature should be considered as an upper limit because of its uncertainty on the low- 
temperature side. Indeed, the finite resolution of the images tends to smear out the sharp 
structure associated with Mott insulator boundaries, leading to an overestimation of the 
actual temperature. To overcome this limit, the spin-gradient thermometry proposed 
in [17] could be employed. 

Summary and concluding remarks 

To summarize, we have shown on various examples of Fermi and Bose gas systems 
how in situ absorption images can provide the grand-canonical equation of state of the 
homogeneous gas. This equation of state is obtained up to a global shift in chemical 
potential and we have given several examples for its determination. The method relies 
on the local density approximation, which is satisfied in many situations, but notable 
exceptions exist such as the case of the ideal Bose gas. The equation of state given by this 
procedure allows direct comparison with many-body theories. While we have illustrated 
here this method on a single-component Bose gas and a spin-balanced Fermi gas, it can 
easily be generalized to multi-component gases. For instance the phase diagram and 
the superfiuid equation of state of spin-imbalanced Fermi gases have been obtained in 
[H [H]. We expect this method to be very useful in the investigation of Bose-Bose, 





knT = 0.09 



1+0.04 
-0.09 
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Bose- Fermi and Fermi-Fermi mixtures. Finally the equation of state of a Bose gas close 
to a Feshbach resonance may reveal thermodynamic signatures of beyond-mean-field 
behavior in Bose-Einstein condensates j48j . 
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Appendix: validity of local density approximation 

Let us now discuss the validity of local density approximation around the superfluid 
transition in our experiment. Along the z axis, the correlation length ^ diverges around 
the transition point z = Zc according to ^ ~ ^f^K-^ ~ Zc)/zc\~'^^ where u = 0.67 is the 
correlation length critical exponent, directly measured in [19], and in agreement with 
u = {2 — a)/ 3. Local density approximation is expected to become inaccurate in the 
region Zc — Sz < z < Zc + 6z where Sz is given by |3T1 [50] : 

6z^^{z, + 6z), I.e. 6z ^ ZcikpZcy^/^^+^K 

Zc is on the order of the cloud size along z, and is much larger that kp^ which is 
on the order of the inter-particle distance. Given the parameters of our experiments, 
(kpZc)"^^^^'^'^'^ ~ 1% and the size 5z where local density approximation is invalid 
is very small. Given the noise of our data (a few %), the deviation from local 
density approximation is thus negligible. Investigating the critical behavior at the 
superfluid transition, such as measuring the critical exponent a, would be an interesting 
development for this method, as proposed in [50] . 
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